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1 INTRODUCTION 



The study of the distribution of matter on large scales, 
based on redshift surveys of galaxies, provides impor- 
tant constraints on models of cosmic structure forma- 
tion. It has been established that galaxy clusters are not 
randomly distributed but tend to aggregate in larger 
systems, the so called superclusters (cf. Bahcall 1988 
and references therein). Superclusters are the largest, 
isolated and dynamically unrelaxed, due to their size, 
objects in the large scale distribution of matter and thus 
they are ideal probes of the initial conditions that gave 
rise to the distribution of matter on such large-scales. 
One's hope is that they can be used for cosmological 
studies in order to test theories of structure formation 
(cf. West 1989). 

Many authors have claimed that the large scale 
clustering pattern of galaxies is characterized by a 
filamentary and sheet-like distribution (cf. Zeldovich, 
Einasto, & Shandarin 1982; Broadhurst et al. 1990; de 
Lapparent, Geller & Huchra 1991). West (1989), Plio- 
nis, Valdarnini & Jing (1992) and Jaaniste et al (1997) 
studied the morphological properties (shape, size and 
orientation) of galaxy superclusters using Abell/ACO 
clusters and found, based on the simple principal axes 
method, that the vast majority of the superclusters are 
flattened triaxial objects, while Plionis et al. (1992) 
found a predominance of prolate (filament-like) super- 
clusters. Further support to the filamentary case was 
presented recently by Sathyaprakash et al. (1998a), 
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(1998b) analysing the shapes of the overdense regions 
using observational (IRAS 1.2Jy) and N-body simula- 
tion data. 

Studies based on the traditional indicators of clus- 
tering, like the two point correlation function, do not 
deepen our knowledge regarding the morphology and 
physics of structures on large scales. Therefore, many 
authors have been using different approaches, based 
on the geometrical properties of the large scale struc- 
ture, in order to investigate such features. Indeed differ- 
ent methods like minimal spanning trees, shape statis- 
tics (cf. Sahni & Coles 1995 and references therein), 
genus-percolation statistics (Zeldovich et al. 1982; Gott, 
Mellot & Dickinson 1986; Coles & Plionis 1991; Yess 
& Shandarin 1996; Sahni, Sathyaprakash & Shandarin 
1997) and Minkowski functionals (cf. Mecke et al. 1994; 
Kerscher et al. 1997) have been used in order to de- 
scribe the global geometrical and topological properties 
of the matter distribution utilising angular and redshift 
surveys of galaxies as well as N-body simulations. 

In this paper we use the recently completed PSCz- 
IRAS redshift survey in order: 



(i) to investigate whether we can reliably identify 
superclusters and measure their shapes in flux-limited 
galaxy samples, 

(ii) to measure the shape and size distribution of the 
PSCz superclusters and 

(iii) to investigate whether these distributions can be 
used as a cosmological probe. 
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The plan of this paper is the following: In section 2 
we describe the PSCz data that we used, in Section 
3 we present the supercluster identification procedure, 
tests of systematic biases that enter in such a procedure 
and the detected PSCz superclusters while in section 4 
we present the supercluster shape determination proce- 
dure, systematic effects that affect their shapes and the 
PSCz supercluster shapes. In section 5 we compare the 
PSCz results with the corresponding ones of three cos- 
mological models and finally in section 6 we draw our 
conclusions. 



2 THE PSCZ GALAXY SAMPLE 

We use in our analysis the recently completed IRAS 
flux-limited 60-/im redshift survey (PSCz) which is de- 
scribed in Saunders et al. (2000) . It is based on the IRAS 
Point Source Catalogue and contains ~ 15500 galaxies 
with flux Su m > 0.6 Jy covering the 84% of the sky. The 
subsample we use, defined by a limiting galaxy distance 
of 240 h^ 1 Mpc, contains ~ 11823 galaxies. 



2.1 PSCz selection function 

Due to the fact that the PSCz catalogue is a flux-limited 
sample there is the well-known degradation of sampling 
as a function of distance from the observer (codified by 
the so called selection function). Unless one corrects for 
this effect the derived galaxy density field will have little 
to do with the true one. This is usually done by weight- 
ing each galaxy by the inverse selection function, (r), 
assuming that the unobserved galaxies are spatially cor- 
related with the observed ones. Note that according to 
the above the selection function is defined as the frac- 
tion of the galaxy number density that is observed above 
the flux limit at some distance r. Therefore 



r=-p-(l-(l + z- Sz)- 1/2 ) (1 + z- 8z) 

where H a = 100 h km s _1 Mpc -1 and Sz is a non- linear 
term to correct the redshifts for the galaxy peculiar ve- 
locities: 

Sz= i(u(r) -u(0)) -f (3) 

with u(0) the peculiar velocity of the Local Group and 
u(r) the galaxy peculiar velocity at position r, which 
can be found in linear theory by: 



3/2 
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where L m in(r) = 4irr vSum is the luminosity of a 
source at distance r corresponding to the flux limit Su m , 
v = 60-/u,m and (p g ) is the mean galaxy number density, 
given by integrating the luminosity function over the 
whole luminosity range, with L m i n = 10 8 h 2 L & (since 
lower luminosity galaxies are not represented well in the 
available samples; cf. Rowan- Robinson et al 2000), and 
10 13 h 2 L @ . Obviously, <p{r) is a decreasing 
function of distance because a smaller fraction of the 
luminosity function falls above the flux limit at greater 
distances. In this work we used a luminosity function 
of the form assumed by Saunders et al. (1990) with 
L« = 10 8 ' 45 h 2 Lq, a = 0.711, a = 1.09 and C = 0.0308 
(cf. Rowan- Robinson et al. 2000). 



2.2 PSCz galaxy distances from redshifts 

It is well known that the distribution of galaxies in red- 
shift space is a distorted representation of that in real 
comoving space (Kaiser 1987). The redshift distance of 
each galaxy can by found from: 



u(r) = 



H /3 
4tt 



S*(r)-. 



:d 3 r 



(4) 



where (3 — Q°' 6 & _1 , assuming linear biasing between the 
background matter and the galaxy fluctuation field 8 g . 
Note that heliocentric redshifts are first transformed to 
the Local Group frame using cz ~ czq +300 sin(Z) cos(b) 
and then the distance to each galaxy is usually found 
by iteratively solving the above set of equations. Such a 
procedure was used by Branchini et al (1999) to recover 
the true density field of IRAS-PSCz galaxies in order 
to study the galaxy density and velocity fields. In our 
work we use their results (for (3 = 1) which were kindly 
provided to us by Dr. Branchini. 



2.3 Incomplete Sky Coverage 

In order to produce a continuous and whole sky density 
field it is essential to treat the 16% of the sky which 
is devoided of data (galactic plane, high cirrus emission 
areas and unobserved regions). We use the PSCz data of 
Branchini et al. (1999), in which they followed the Yahil 
et al. (1991) method to fill the galactic plane region with 
synthetic objects that reproduce the mean density of 
galaxies in the regions residing nearby. At higher galac- 
tic latitudes they fill the masked regions again with syn- 
thetic objects so that they reproduce the overall mean 
number density of PSCz galaxies. 



3 IDENTIFYING SUPERCLUSTERS IN 
THE SMOOTH PSCZ DENSITY FIELD 

3.1 Smoothing Procedure 

For the purpose of this study we need to derive from the 
discrete distribution of PSCz galaxies a smooth contin- 
uous density field. To this end we use a Gaussian kernel 
on a N 3 grid: 



W(| Xj -x g |) = (2tvR 2 s 



-3/2 



exp - 



2RL 



.(5) 



The smoothed density, at the grid-cell positions x g , is 
then: 

£iPte)W(l*i-x«|) 

J MI x j - x g |)d 3 a; 

where the sum is over the distribution of galaxies with 
positions x 3 . 

In our study we will use two smoothing radii, 
namely i? sm = 5 h^ 1 Mpc and 10 h^ 1 Mpc, to probe 
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essentially different supercluster sizes. We will extend 
our analysis out to r max , where r max — 150 and 240 

Mpc respectively for the two smoothing radii while 
the size of each cell is set equal to 7? sm . Therefore we 
use a grid of size of 60 3 and 48 3 respectively. Note also 
that |xj — x g | < 3i?sm and thus the integral in the de- 
nominator of eq. (^|) has a value ~ 0.97. 

In order to construct our supercluster candidates 
we select all cells with an overdensity above a cho- 
sen threshold and join together those having common 
boundaries. Due to the fact that superclusters should 
be identified as very high density peaks, the critical 
value of the overdensity threshold is defined directly 
from the probability density function as its the 97%- 
quantile value (using lower values we tend to percolate 
structures through the whole volume). 

3.2 Correcting Biases Introduced by 
Smoothing 

Due to the drop of the selection function as a func- 
tion of distance, shot-noise effects increase with dis- 
tance. Although Gaztanaga & Yokoyama (1993) have 
shown that the smoothing process itself considerably 
suppresses such shot-noise effects, the coupling between 
the selection function and the constant radius smooth- 
ing will result in a distorted smoothed density distribu- 
tion, especially at large distances. We expect that Gaus- 
sian spheres, centered on distant cells, will overestimate 
the true density in regions where galaxies are detected 
(due to the large <j>{r) weighting), while in underdense 
regions they will underestimate the true density. 

The decrease of resolution, due to the drop of the 
selection function at large distances, can be dealt by in- 
creasing the smoothing radius as a function of distance. 
Using the mean inter-galaxy separation as the Gaussian 
smoothing radius is usually a good rule-of-thumb crite- 
rion to suppress such effects in velocity field studies (cf. 
Branchini et al. 1999). However, it cannot be used to 
identify and study supercluster shapes or in studies of 
the pdf and its moments, in which it is necessary to have 
a constant resolution element through out the studied 
volume. 

3.2.1 Quantification of the effect 

We quantify this effect by using N-body simulations for 
which we have the overall 3D "galaxy" distribution as 
well as the corresponding PSCz look-alike distributions 
(a detail study of the effect and of the devised correction 
procedure on the corresponding pdf statistics will be pre- 
sented in Plionis, Basilakos & Rowan- Robinson 2001). 
Here we just note that we use a CHDM (f2hot=0.3 and 
Sib = 0.1) simulation based on an optimize version of 
the Truncated Zeldovich Approximation, described in 
detail in Borgani et al. (1995; 1997). There are 256 3 
grid points and as many particles, while the simulation 
box has L — 480 h~ l Mpc and periodic boundaries. The 
choice of the CHDM model was based solely on the fact 
that the corresponding distribution of particles that we 
tag as 'galaxies' have a correlation length that resembles 
that of the IRAS galaxies (r ~ 4 h~ l Mpc). 




r (h _1 Mpc) 

Figure 1. A(r) as a function of distance for the simulation 
smooth density field with i? S m = 10 ft. -1 Mpc. Filled circles 
represent fluctuation values for p > (p g ) while open circles 
for p < (pg). The continuous lines represents the best poly- 
nomial fit to A(r) and the broken line represents the (A(r)) 
values. 

The main question that we want to answer now 
is: How the 'PSCz'-based smoothed density field relates 
to the corresponding true underlying 3D density field? 
We answer this question by comparing the smoothed 
3D distribution of the parent population of simulated 
"IRAS" galaxies with that of the corresponding PSCz 
look-alikes, generated from this distribution. We derive 
the statistical deviation of the smoothed 'PSCz'-based 
distribution from the corresponding 3D one as a func- 
tion of distance. Therefore, we evaluate the following 
quantity: 



If the selection process does not affect the 'PSCz'-based 
density field then A(r) in the average should be negli- 
gible for all distances. Indeed we find that (A(r)) ~ 
for all distances (broken line in figure 1 & 2). However, 
if we split the densities in two subsamples: one higher 
and one lower than (p g ), then we find a strong distant 
dependent effect (see figures 1 & 2), where the coupling 
between the selection function and the constant radius 
smoothing results in an underestimation and overesti- 
mation of the low and high density regions respectively, 
especially at large distances. 

We can attempt to fit A(r) using a high-degree 
polynomial and then correct the raw smoothed PSCz 
densities by weighting them with the factor W(r) = 
(1 + Aflt(r)) -1 . We have however to note that any such 
multiplicative correction procedure cannot correct (ar- 
tificial) zero densities. To this end we use a six degree 
polynomial and a x 2 minimization procedure to test the 
goodness-of-fit. The fit, presented in figure 1 & 2 as the 
continuous lines, has a reduced x 2 value of < 1 indi- 
cating that the fit is good although the scatter is prob- 
ably underestimated. In Plionis, Basilakos & Rowan- 
Robinson (2001) we will present details on the correc- 
tion procedure and the results in recovering the under- 
lying pdf and its moments. Here we test our procedure in 
recovering the true number, size, multiplicity function 
and shapes of overdense regions (see below). 
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Figure 2. Same as for figure 1 but for the R sm = 5 h 1 

Mpc case. 




Figure 4. Supercluster multiplicity function for the 3D 
(points), raw-'PSCz' (broken line) and corrected-'PSCz' 
(continuous line) cases. 
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Figure 3. Supercluster number density in equal volume 
shells as a function of distance for the 3D (open circles), 
raw-'PSCz' (crosses) and corrected-'PSCz' (squares) cases. 



3.2.2 Robustness of Supercluster Identification 

The above mentioned systematic effect should also af- 
fect the connectedness of overdense regions and thus the 
detected supercluster number, size and shape. To this 
end we compare the number density of superclusters, 
found by our method in the smoothed 3D simulation 
distribution with the average found in 10 'PSCz' real- 
izations. As discussed before, we select cells that have 
an overdensity value above a specific quantile of the pdf. 
However, since we are interested in comparing the su- 
perclusters found in flux-limited samples with the true 
underlying superclusters we have to make sure that we 
are using equivalent thresholds. Such a choice is easy in 
the simulation test case, in which we know a priori the 
true 3D pdf. We have found that in general a specific 
quantile of the raw 'PSCz' look-alike pdf is significantly 
larger than the corresponding 3D pdf an effect that our 
density correction procedure restores to a good degree, 
especially in the pdf high density tail (Plionis, Basilakos 
& Rowan-Robinson 2001; in preparation). Therefore, we 
use as our density threshold a specific pdf quantile. 

In figure 3 we compare the supercluster number 
density in equal volume shells (5V ~ 1.16 x 10 7 ft" 3 
Mpc 3 for the R 3m = 10 ft" 1 Mpc case and SV ~ 
2.83 x 10 6 ft" 3 Mpc 3 for the R sm = 5 ft" 1 Mpc case) 



as a function of distance for the 3D, raw-'PSCz' and 
corrected-'PSCz' distributions. For the R sm = 10 ft" 1 
Mpc case we find that the mean supercluster number 
density per equal volume shell is 17±13, 10.9±3 and 
9.3±3.3 for the 'raw', 'corrected' and true 3D density 
fields, while the corresponding densities for the R am = 5 
ft" 1 Mpc case are 30±13, 13±4 and 14±3. 

It is evident that for the 'raw-PSCz' case the su- 
percluster number density increases artificially with dis- 
tance; at large distances the galaxy sampling rate is low 
and the smoothing procedure is unable to produce a 
continuous density field, producing a large number of 
artificially disjoint regions. This effect is efficiently sup- 
pressed once we correct the smoothed density distri- 
bution according to our correction procedure for both 
smoothing scales used. 

In figure 4 we present the multiplicity function of 
the identified superclusters. Again we see that our cor- 
rection procedure is effective in recovering the underly- 
ing multiplicity function. 

We conclude that the connectedness of overdense 
regions in flux limited samples, with a selection func- 
tion similar to that of IRAS 0.6 Jy galaxies, can be 
recovered once we correct the 'raw' smoothed density 
field according to our procedure. We have verified, how- 
ever, that in order to restore such connectedness we need 
to correct the smoothed density field limiting ourselves 
within ~ 230 ft -1 Mpc for the R sm > 10 ft" 1 Mpc case 
and within ~ 150 ft" 1 Mpc for the R am > 5 ft" 1 Mpc 



3.3 Space Density and Multiplicity of PSCz 
Superclusters 

In figure 5 we present the space density of PSCz super- 
clusters in equal volume shells, having size as that in the 
previous section, before and after correcting the density 
field according to our procedure. We have also verified 
that variations of 8V do not alter our results. For the 
R sm = 10 ft" 1 Mpc case, in which we are revealing 
relatively large structures, we find in the uncorrected 
density field that the first shell (< 130 ft" 1 Mpc), con- 
tains no supercluster, which evidently is false since it is 
well known that within this distance there are several 
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Figure 5. Mean number density of PSCz superclusters in 
equal volume shells, based on the corrected (solid circles) 
and uncorrected (open circles) smooth density fields. 




Figure 6. Multiplicity function of PSCz superclusters for 
the corrected (continuous line) and uncorrected (broken line) 
smooth density fields. 



superclusters like the Great Attractor, Perseus-Pisces, 
Coma etc. We also observe the systematic increase of 
supercluster space density with distance in both R 3m 
cases. However, once we correct the density field ac- 
cording to our procedure, discussed in 3.2, we recover a 
constant space density of superclusters in the different 
equal volume shells. 

In figure 6, we present the PSCz supercluster mul- 
tiplicity function, for the uncorrected (dashed line) and 
corrected density fields (continuous line), which will be 
used in the last section as a cosmological probe. In table 
1 and 2 we present for both smoothing scales, the de- 
tected PSCz superclusters, having a volume above the 
indicated threshold, as well as their shape parameters 
(see section 4). We have identified most of the known 
superclusters but also quite a few new ones. Many other 
known superclusters, not shown in the tables, have been 
found but fall below the specific volume threshold. 



4 3D SUPERCLUSTER SHAPES 

4.1 Shape Determination Method 

The supercluster shapes are defined by fitting ellipsoids 
to the data. Shapes are estimated for those "superclus- 
ters" that consist of 8 or more cells. To this end we use 
the moments of inertia method (cf. Carter & Metcalfe 



1980; Plionis, Barrow & Frenk 1991): In = w i x i, 

hi = W iVil hi = -^21 = J2 W i X iVi> ^3 = hi = 

WiXiZi, I23 = I32 = Wii/iZi, with Wi the statis- 
tical weight of each cell. This is defined as the density 
fluctuation: 

„., _ Pi(Xgr) - (fig) 



(P 3 ) 



(8) 



where (p s ) is the mean galaxy number density. Note 
that because of its symmetry we diagonalize the inertia 
tensor 

det(Jy - A 2 Af 3 ) = (M 3 is 3 x 3 unit matrix) , (9) 

obtaining the eigenvalues I\, h, h from which we de- 
fine the shape of the configuration since the eigenvalues 
are directly related to the three principal axes of the 
fitted ellipsoid. The corresponding eigenvectors provide 
the direction of the principal axes. 



4.2 Geometrical Shape Statistics 

To identify the characteristic morphological features of 
superclusters and to investigate the geometrical prop- 
erties of the large-scale structure, we will use the dif- 
ferential geometry approach, introduced by Sahni et al. 
(1998). 

We remind the reader the basic elements of differ- 
ential geometry theory (cf. Lipschutz 1969). Assuming 
a local coordinate system on a surface with the inte- 
grating class C 2 , it is well known that the geometrical 
features of the surface are described by the first and 
second fundamental forms: 



I = dr dr = Ed6 2 + 2Fd6d(f> + Gdcf> 2 



II 



-dr ■ dn = Ld6 z 



2Md6d4> + Nd(f> 2 



(10) 
(11) 



where E = rg ■ rg, G — r^ ■ r^ F = r e ■ r^, L = rgg ■ n, 
G = r^^ • n and M = r^ ■ n H (n is the perpendicular 
vector n = reXT * ) Thus, the surface area, integrated 
mean curvature and genus can be written: 



S = 



c 



\J EG - F 2 d8d(j} 



ki + fej 



dS 



47T 

where 



fci + ki — 



2 

fcifcadS . 

EN + GL- 2FM 
EG — F 2 



and 



LN — M 
EG — F 2 



(12) 
(13) 
(14) 

(15) 

(16) 



The genus parameter characterises the surface topology; 
multiply connected surfaces have Q > while simply 
connected ones have Q < 0. 



~rr — cllld TCrvn 

a a. 



() a. - 



© 0000 RAS, MNRAS 000, 000-000 



6 Basilakos, et al. 



Table 1. Detected superclustcrs in the R, 
distance < 210 h^ 1 Mpc. 

Volume (10 4 h,-' 3 Mpc 3 ) (d) 



41 


146. 


70 


177. 


36 


195. 


30 


148. 


38 


201. 


27 


152. 


103 


171. 


51 


136. 


33 


110. 


40 


159. 



= 10 ft" 1 Mpc density field with V 







Z sup 


K X 


K 2 


40. 


-140. 


-13. 


.095 


.113 


-85. 


-102. 


-118. 


.104 


.076 


10. 


-101. 


-167. 


.067 


.126 


-107. 


-96. 


35. 


.084 


.095 


-187. 


-75. 


-10. 


.012 


.016 


104. 


-55. 


-95. 


.018 


.020 


-157. 


58. 


38. 


.149 


.098 


-117. 


70. 


-4. 


.049 


.083 


-32. 


75. 


74. 


.092 


.065 


-55. 


95. 


-115. 


.043 


.066 



2.5 x 10 4 h~ 3 Mpc 3 (N cell > 25) and 



K 1 /K 2 Name 
.839 Pisces 
1.368 
.534 
.886 
.745 
.888 

1.510 Far Shapley 
.585 Near Shapley 

1.415 Hercules 
.645 Near Horologium 



Table 2. Detected superclusters in the R sm = 5 h 1 Mpc density field with V > 3.5 X 10 3 h 3 Mpc 3 (N cen > 28) and distance 
< 140 h- 1 Mpc. 



;io 3 h- 3 


Mpc 3 ) 


(d) 




Ysup 


Z S up 


K X 


K 2 


Ki/K 2 


Name 


4.1 




134. 


46. 


-125. 


0. 


.054 


.118 


.458 


Pisces 


4.1 




137. 


-94. 


-96. 


31. 


.044 


.058 


.764 




6.0 




108. 


59. 


-78. 


48. 


.057 


.033 


1.697 




6.0 




105. 


32. 


-58. 


-81. 


.112 


.215 


.522 




3.5 




123. 


52. 


-46. 


102. 


.111 


.048 


2.328 




4.9 




89. 


72. 


-47. 


-24. 


.054 


.082 


.655 




5.0 




128. 


26. 


-36. 


-120. 


.060 


.074 


.805 


Lcpus 


5.4 




138. 


95. 


-46. 


-90. 


.049 


.054 


.908 




7.0 




49. 


-39. 


-24. 


18. 


.078 


.113 


.692 


Pavo-Indus 


7.0 




49. 


46. 


-14. 


-9. 


.116 


.053 


2.205 


Perseus- Pegasus 


15.1 




44. 


-36. 


19. 


-15. 


.080 


.055 


1.460 


Great Attractor 


6.3 




137. 


-118. 


49. 


50. 


.029 


.037 


.779 




11.8 




132. 


-114. 


67. 


-1. 


.071 


.158 


.445 


Near Shapley 


3.9 




107. 


87. 


52. 


-32. 


.000 


.000 






5.5 




110. 


-43. 


77. 


65. 


.055 


.138 


.400 


Hercules 


9.0 




119. 


-25. 


75. 


89. 


.090 


.196 


.459 




3.5 




78. 


3. 


77. 


14. 


.002 


.002 


.851 


Coma/A1367 


7.3 




124. 


-14. 


107. 


61. 


.063 


.106 


.591 




6.8 




126. 


-44. 


108. 


-47. 


.033 


.063 


.517 




3.6 




137. 


4. 


122. 


-63. 


.045 


.108 


.418 





Furthermore, we define the 3 dimensional shape of 
structures by fitting the best triaxial ellipsoid (M — 0), 

r(8, 4>) — Jisin#cos0i + ^sin^sin^j + Iscos9k , 

to the different isodensity contours having volume V = 
^-hhh and < 4> < 2tt, < 9 < n. 

Here we present the basic steps of the shape statis- 
tics that we will use, following the notation of Sahni et 
al. (1998). They introduced a set of three shapefinders 
Tii = VS~\ Ti 2 = SC- 1 and Ti 3 = C, having dimen- 
sions of length. In the same framework they defined two 
dimensional shapefinders K = (Ki,K 2 ), where 



Ki = 



and 



K 2 = 



TL 2 — Tii 
Ti 2 + Tii 



TLz — Ti 2 



(17) 



(18) 



Tiz + Ti 2 ' 

normalized to give Tii — R (Ki, 2 — 0) for a sphere 
of radius 7?. Therefore based on these shapefinders we 
can characterise the morphology of cosmic structures 
according to the following categories: 

• Pancakes for Ki/K 2 > 1 

• Filaments for Ki/K 2 < 1 



• Ribbons for (Ki,K 2 ) — (a, a) with a < 1 and thus 
Ki/K 2 ~ 1. 

• Sphere for h ~ I 2 ~ I 3 and thus (Ki,K 2 ) ~ (0, 0). 

• Ideal filament (1-d objects) for Tii ~ Ti 2 << Tiz 
and thus (K U K 2 ) ~ (0, 1). 

• Ideal pancake (2-d object) for Tii << Ti 2 ~ Tiz 
and thus (Ki,K 2 ) ~ (1,0). 

Finally, for the quasi-spherical objects the Ki <2 are 
small and therefore their ratio (Ki/K 2 ) measures the 
deviation from spherical shapes. 



4.3 Robustness of Supercluster Shape 
Determination 

Here we investigate how do the systematic effects, dis- 
cussed in section 3.2, affect the determination of super- 
cluster shapes. To this end we compare the shapes of the 
superclusters, found by our method in the smoothed 3D 
simulation distribution with those found in 10 'PSCz' 
realizations. 

Performing many tests using different smoothing 
radii, we have found that the ratio of the dimensionless 
supercluster shapefinders ( "shape spectrum" ) is also af- 
fected from the bias discussed above. Note that we ex- 
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Figure 7. Shape Spectrum fluctuations between the 3D 
and raw-'PSCz' case (crosses) and between the 3D and the 
corrected-'PSCz' case (squares). 
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Figure 8. The corrected (hatched region) and uncorrected 
(dashed line) shape spectrum of PSCz superclusters for the 
two smoothing radii. 



elude the ideal spherical surfaces, Ii ~ I2 — I3, be- 
cause their ratio K\ /K2 tends to infinity. In figure 9 we 
present the fluctuations of the shape spectrum between 
the 3D and raw-'PSCz' case and between the 3D and the 
corrected-'PSCz' case. Again it is evident that our cor- 
rection procedure suppresses significantly the system- 
atic effects introduced by the coupling between the se- 
lection function and the constant radius smoothing for 
both smoothing radii. 

We conclude that in order to correctly identify su- 
perclusters and their connectedness, as well as their 
shapes from flux limited samples, it is necessary to take 
into account the bias of the smoothed density field, 
discussed in section 3.2, and correct accordingly the 
smoothed galaxy distribution. 

4.4 The PSCz Supercluster Shape Spectrum 

In figure 10 we present the uncorrected and corrected 
"shape spectrum" for both smoothing radii. It is obvi- 
ous that the dominant feature of the PSCz supercluster 
shapes is filament ariness; ie., K1/K2 < 1 (in agreement 
with previous studies; see Sathyaprakash et al 1998). 

Furthermore, although the fraction of filamentary 
superclusters (K1/K2 < 1) remain roughly the same be- 
fore and after the correction, there is a slight change of 
the corrected 'shape-spectrum' towards less extreme fil- 
amentary superclusters. Regarding extreme shaped su- 
perclusters (with N ce u > 8), we have found, in the 
Rsm — 5 k" 1 Mpc smoothed field, 8 very filamentary su- 
perclusters with K1/K2 < 0.45 (among which the Near 
Shapley and Hercules superclusters), 6 spherical and 
only 2 extreme pancake-like structures with K\ / K2 > 2 
(one of which is the Perseus- Pegasus supercluster; see 
table 2). In the R sm = 10 h -1 Mpc smoothed field we 
find only one such extreme filamentary supercluster. 



5 COMPARISON WITH COSMOLOGICAL 
MODELS 

We use mock PSCz catalogues, generated from three 
large cosmological N-body simulations, in order to in- 
vestigate whether supercluster properties (shapes and 



numbers) can discriminate between models. The cos- 
mological simulations are based on an AP 3 M code and 
they are described in Cole et al. (1998). We consider 
three different cold dark matter models: (1) a flat low- 
density CDM model with tt m = 1 - Q A = 0.3, (2) a 
critical density universe fl m = 1 with power spectrum 
F = 0.25 (rCDM) and (3) a critical density universe 
n m = 1 with power spectrum T = 0.5 (SCDM). The 
first two models are normalized by the observed cluster 
abundance at zero redshift; having fluctuation ampli- 
tude in 8 /i _1 Mpc scale of ct 8 = 0.55fi~°' 6 (Eke, Cole, 
& Frenk 1996) , while the third is COBE normalized with 
(78 = 1.35. 

For each cosmological model we average results over 
10 nearly independent mock PSCz catalogues extend- 
ing out to a radius of 170 ft _1 Mpc (cf. Branchini et al. 
1999). Good care has been taken to center the cata- 
logues to suitable LG-like observers (having similar to 
the observed Local Group velocity, shear and overden- 
sity). 

We analyse the mock PSCz catalogues density fields 
(with Rsm — 5 ft -1 Mpc) in the same fashion as that of 
the observed PSCz catalogue and we compare the out- 
come of this procedure in figure 11, where we plot the 
detected supercluster shape-spectrum and multiplicity 
function for all three models and PSCz data. Evidently 
the SCDM model performs worst in the comparison. In 
order to quantify the differences between models and 
data we perform a standard \ 2 test and we present the 
V >x 2 results in table 3. The only model that is excluded 
by the data at a relatively high significance level is the 
SCDM model, while the ACDM model best reproduces 
the observed supercluster shape-spectrum and multi- 
plicity function. The rCDM is marginally excluded by 
the multiplicity function comparison but only at a ~ 
94% level. 

Comparing the models among themselves we see 
that the shape-spectrum is insensitive to the different 
cosmologies, probably because supercluster shapes re- 
flect the Gaussian nature of the initial conditions which 
are common to all models. However, the supercluster 
multiplicity function is a strong discriminant between 
the models and therefore our results do show a prefer- 
ence of the ACDM model over the Q m = 1 models that 
we have studied. 
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Table 3. x probabilities (J } > 2) of consistency between 
data and models as well as between the models themselves. 



Comparison Pair 
PSCz - ACDM 
PSCz - SCDM 
PSCz - tCDM 

ACDM - SCDM 
ACDM - rCDM 
tCDM - SCDM 



Shape-Spectrum 
0.148 
0.019 
0.146 

0.71 
0.94 
0.07 



Multiplicity fun. 
0.135 
0.019 
0.06 

5.2xl0" 3 
6.6xl0" 4 
2xl0~ 8 
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Figure 9. PSCz Superclusters comparison with Cosmolog- 
ical models: Left panel: Multiplicity function, right panel: 
Shape spectrum. The PSCz results are represented by the 
hatched region. 



6 CONCLUSIONS 

We have studied the properties of superclusters detected 
in the smoothed PSCz galaxy density field as connected 
regions above some overdensity threshold. We have in- 
vestigated, using simulations, the biases that enter in 
the detection and shape determination of superclusters 
which are identified in smoothed density fields of flux- 
limited galaxy samples. 

We have devised a statistical approach to correct 
for such biases which we have verified that indeed it 
recovers the underlying 3D properties of superclusters. 
To determine supercluster shapes we use the differen- 
tial geometry approach of Sahni et al. (1998) and find 
that the dominant oupcrclustcr morphological feature ia 
filamentarinasB , 



Finally, we have compared our PSCz supercluster 
results with the corresponding ones generated from the 
analysis of three cosmological models (SCDM, r-CDM 
and ACDM) and we find that the model that best re- 
produce the observational results is the ACDM model 
(«A = 0.7). 
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